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Abstract 

We prove that a stochastic process of pure coagulation has at any time t > a time 
dependent Gibbs distribution if and only if the rates of single coagulations are of 

the form = if(j) + jf{i)i where / is an arbitrary nonnegative function on the 

set of positive integers. We also obtain a recurrence relation for weights of these Gibbs 
distributions that allow to derive the general form of the solution and the explicit solutions 
in three particular cases of the function /. For the three corresponding models, we study 
the probability of coagulation into one giant cluster, by time t > 0. 
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1 Process of pure coagulation: formulation 
of the model. Objective and the context 
of the paper. 

We consider a standard model of stochastic coagulation (see e.g. [7]), viewed as a time 
continuous Markov chain on the finite set Qn of partitions rj of a given integer iV : 

N 

tt N = {q = (ni, . . . , n N ) : ^ kn k = N}. 

k=l 

In the context considered, N is the total number (=total mass) of identical particles 
(molecules, planets, animals etc) partitioned into clusters of different sizes, so that n k is 
the number of clusters of size k in a partition 77 G Ojv- Infinitesimal in time transitions 
are coagulations of any two clusters of sizes i and j into one cluster of size i + j, resulting 
in the state transition 

77 ->■ 

Here 77 g £In codes the state that is obtained from a state 77 G fijv, with > 0, nj > 0, 
by a coagulation of any two specific clusters of sizes % and j. In the sequel such coagulations 
are called single and their rates are denoted The following assumptions on the 

rates ip(i,j) and on the induced rates ^(77 — > 77^) of state transitions, describe the class 
of coagulation processes (CP's) considered: 

• The function ip(i,j) is nonnegative, symmetric in i,j and is not dependent on iV; 

• The rate ^(77 — > 77^) is equal to the sum of rates of all single coagulations 
of rii > groups of size % with rij > groups of size j each one, so that 

K(rj -> 77 (ij) ) = mnrffaj), i j, 2<i + j<N, 
K( V ^ V ^) = Ul{jlt ~ l) 2<2i<N. 

An interpretation of the above expression for the coagulation kernel K in terms of 
the kinetics of droplets of different masses is given in [20J. CP's with rates of state 
transitions of the above form are naturally called mean-field models, meaning that at any 
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state rj G Qn, any cluster can coagulate with any other. Hence, given N, the distribution 
of a CP = CP(N) at any time t > is uniquely specified by the initial distribution on 
the set Qn and the rates ip{i,j). 

The history of CP's goes back to 1918 when Smoluchowski formulated a deterministic 
version of the model of pure coagulation of molecules in chemical kinetics in his seminal 
paper. Because of the ever growing field of applications and a rich probabilistic context, 
the study of a variety of versions of the model continues to be a hot topic in the theory 
of stochastic processes. Marcus [20] was apparently the first to formulate the stochastic 
version of a CP. A particular case of a pure coagulation, when j) = a(i+j)+b, a, b > 
is known in the literature as the Marcus- Lushnikov stochastic process, while in [T] as well 
as in some other papers, the name is given to all stochastic CP's with rates of single 
coagulations of the form N~ x i])(i) j), with an arbitrary if)(i,j). It is important to point 
out that, in contrast to the latter Marcus- Lushnikov process, the basic assumption of our 
setting is that the rates of single coagulations do not depend on N. The equilibria of 
some reversible models with rates of coagulation and fragmentation depending on N were 
studied in [13], [H]. 

Let Xff(t), t > denote a CP(N) starting from an initial distribution p on Qpf, The 
objective of our study is the probability distribution (=transition probability) p(r], p; t) of 
the process, which is 

F(X$(t) = v ), v eQ N , t>0. (l.i) 

In the sequel we refer in more details to the literature related to the aforementioned 
objective. 

We describe now the organization of the paper. In Section 2 we formulate our main 
result which is the characterization of CP's possessing a probability distribution (II. ip of 
Gibbsian type at any time t > 0. As corollaries of our theorem we derive the general 
form of weights of the aforementioned Gibbs distributions and the explicit expressions 
for the weights for three particular models of CP's. We also analyze the behaviour in 
time of some important functionals of the models. In the final section, Section 3, we 
explain the name "Gibbsian" given to the class of distributions considered, describe the 
linkage to coagulation-fragmentation processes on set partitions and indicate the relation 
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of Gibbsian distributions to the theory of random combinatorial structures. 

2 Main result 

Recall that X£ (t), t > denotes a CP(N) starting from an initial distribution p on Qjy. 
Our goal is to identify CP's with a probability distribution ( 11. ip of the form 

*, P ;') = c„(i)n^ 

k=l Uk ' 

7] = (m, . . . ,n N ) G Q N , t>0. (2.2) 

Here Cjv(t), t > is time-dependent partition function and a^jy^) > 0, t > are 
time-and iV-dependent weights. The distributions in the right hand side of (12.21) are 
called Gibbsian (=canonical Gibbsian). We note that by the above definition the initial 
distribution p is Gibbsian. By virtue of the mass conservation law YljLi 3 n j = tilting 
transformations of the weights in (12.21) with an arbitrary function r/v(i) > 0, t > : 

a k ,N(t) = (r N (t)) k a k ,N(t), k = l,...,N, t>0 

and the induced transformation C^(t) = (riv(t))~ Ar Civ(t) of the partition function, result 
in different representations of the generic Gibbs distribution. In view of this fact we 
assume in the rest of the paper that the weights ak^if) in (12. 2p are such that the partition 
function does not depend on t > 0, i.e. CV(t) = CV, t > 0. The latter assumption appears 
to be of great help for our discussion below. 
Our main result is the following characterization. 

Theorem. The probability distributions p(rj , p; t) of a CP X^\t), t>0 are of a Gibbsian 
form 112. if and only if the following three conditions are satisfied: 

(i) The initial distributions p on Q N are Gibbsian with arbitrary weights a k<N = 
a fc ,jv(0) > 0, k — 1, . . . ,N, N > 1. 

(ii) The rates of single coagulations are of the form 



Hhj) = if(j)+jf(i), i,j:l<i + j<N, 



(2.3) 
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where f is an arbitrary nonnegative function on the set of positive integers. 
(Hi) The weights a kjN (t) are defined recursively by 

a hN (t) = e-C- 1 "'"^, t > 0, 
Jo 2 

k = 2,...,N, t>0, (2.4) 
where > are constants implied by the initial distribution p in (i). 

Proof. Our plan is to show that the assertions (i) — (Hi) of Theorem are implied 
by the Kolmogorov forward equations. Suppose is the state that is obtained from a 
state rj G Qn, with n^+j > 0, by a fragmentation of some cluster of size i+j>2 into two 
clusters of sizes i and j. Then the equations read as follows: 



dt 



p(v,p;t) = -p(v,p;t)( niU 3 "^(ijj) + ^(M)) + 



l<i<j<N i=l 



l<i<j<N 



/- ^ (rif + l)(nj + 2) 
l^P(V(i,i), P,t) « ( 2 - 5 ) 



»=i 



First, we assume that (12.21) holds. Substituting (12. 2p with weights afc,jv(t), such that the 
partition function does not depend on t > 0, gives 

l<i<j<N i=l 

. / ^ a^t)a^) ( 1)( ^.^ - )+ 

i<fe<iv a *+i.^( t ) + 1 )( n J + X ) 

a-,jy(t) n 2i (n i + l)(n i + 2) 

+ l)(n i + 2) 2 

which can be written as 



• rw ' K ' y l<i,j<N l"WVW l<i<j<N i=l 
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rieQ N , t>0. (2.6) 
In view of ( 12. 2 p and the time independence of the partition function, we proceed as 

E X ai,jv(*) \^ a i7N (t)a j>N (t) / ^ ,/■ ,,p^(ni-l) ,/. .x 

fc=l i,JVV 7 l<i,j<N l +J. Jv v / l<i<j<N i=l 

r] e Q N , t > 0. 

Finally, we have 

E^f* = (L L ^M-^( M H) 

fe=l vvv ; fc=2 ft ' JVV ; l<i,j<N 8=1 

77 G fytf-, * > 0. 

We now rewrite the last equation as 
1 ^ 

-( E ^j^j • tpjiJ)) = ^ A k ,N(t)n k , t > 0, for all 77 = (m, • • • ,n N ) € fW, (2.7) 

l<i,j<N k=l 

where 

^fc,M*) = " o 77\ —F7\ + o^( k > k )> k = l,...,N, (2.8) 

assuming that for k = 1 the sum in the right hand side of the last expression equals zero. 
Applying the mass conservation law, the equation (12.71) conforms to 

\( Yl UiU i ■V'(^i)) = E ninj(jA iiN (t)+iA jtN (t)), 

^<i,3<N l<i,j<N 

for all 7] = (ni, . . . , n N ) e Q N , t > 0, 

from which, using the assumed symmetry of the function we derive the unique 

form of ij) : 

*l)(i,j) = ^-(jA iiN (t)+iA j , N (t)), (m,...,n N ) eQ N , t>0. (2.9) 
Consequently, 

2 

= —iA i>N (t), t>0, i>l, 
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which imposes the following necessary and sufficient conditions on the coefficients A iN (t) : 
iV _1 A )iV (t) =: f(i), i>l does not depend on JV > 1 and t > 0. (2.10) 
By virtue of this condition and (12.91) . we find 

il>{i,i) = 2if(t), 1>(iJ)=ifti)+3fV), (2- 11 ) 

which proves the necessity of part (ii) of the theorem. It is left to demonstrate that the 
weights a^N^t) can be found recursively from (I2.10p and (12.81) . Firstly, we recover ai^t) : 

n-\- +/(i)) = mi t>o, 

ai >N {t) 

ai, N (t) = a 1:N e~ {N ~ mi)t , t > 0. (2.12) 

where aijv(0) = a^jv > is a constant, given by the initial distribution p. For k > 2, 
fl2~TUD takes the form 

V 2a ktN (t) a k ,N{t) ' 

Observing that the convolution term 

M k)N {t) := ^2 a i,N(t)a j>N (t)tJj(iJ)= ^ a i,N(t)a jtN (t)(if(j)+jf(i)), t>0 

i+j=k i+j=k 

does not depend on a^jsfit), we arrive at the following first order differential equation for 

a[ jN (t) = ^f^--(N-k)f(k)a k , N (t), k = 2,...,N. (2.13) 
Solving it we obtain the recurrence relation that conforms to part (Hi) of the theorem: 

a k>N (t) = e -("-*)/W ( y* MMM. e (N - k ^ k >du+a^ N ) , k = 2, . . . , N, t > 0. (2.14) 

Now the validation that, under conditions (i) — (Hi), (12. 2p is the solution of the Kol- 
mogorov equations (I2.5P is obvious. ■ 
Notes. (1) An obvious consequence of the theorem is that a CP with given rates i/i(i,j) 
of form (I2.3p . starting from a non Gibbsian initial distribution, does not possess Gibbsian 
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transition probabilities. In particular, if the aforementioned process starts from a mix- 
ture of Gibbs distributions, then its transition probabilities ( II .ip are the corresponding 
mixtures of Gibbs distributions. 

(2) The differential recurrence (12.41) was originally derived by Lushnikov in |18j . where 
the weights a^A^t), t > are associated with the generating function for the probabilities 
p(r], p;t), t > 0. Buffet and Pule [S] proved the core fact that, given ifj(i, j) = if(j) +jf{i), 
the Kolmogorov equations are solved by Gibbsian distributions p(r], p\ t) with weights sat- 
isfying the differential recurrence (12.41) . Our result strengthens the result of [5] in the 
following four directions, (i) We establish the necessity of the form (12.31) of rates for 
Gibbsian transition probabilities (II. ip ; (ii) The study in [2] is limited to Gibbs distribu- 
tions with partition functions Cat = 1, N > 1, which substantially restricts the class of 
possible initial distributions. For example, under the above restriction, the initial Gibbs 
distribution with constant weights a^jy = 1 is not allowed, since in this case the explicit 
expression for Cjy as a function of N is not known. So, the trick with tilting transfor- 
mation can not be applied, (in) Our theorem is proven for the case of a nonnegative 
(rather than positive) function /, which enables us to treat CP's like those of Becker- 
Doring, which are defined in Corollary (D) below; (iv) Solving in Corolary (A) below the 
recurrence relation (12. 4p . we find a general form of the weights cik,N(t). 

(3) In |5j, p. 1047 it was noted that if N is a multiple of an integer q > 1, then the 
initial distribution p concentrated on the partition of N into y groups of size q each one, 
is Gibbsian, with 

a q>N > 0, a fcjA r = 0, N > k q. 

An important particular case of the above measures is given by q = 1, which in chemistry 
corresponds to what is called total dissipation. In this case the initial distribution p is 
concentrated on the partition (N, 0, . . . , 0) and CV = -w~, where a± jy > is arbitrary, 

a l,JV 

while cik,N = 0, k — 2, . . . , N. A closer look on the above example in [5] shows that a class 
of Gibbs initial distributions encompasses a variety (but not all) of measures concentrated 
on single partitions of N. We give an example of one of such measures. Let N — I + m, 
where < I < m and iV is not divisible by I. This secures that rj* = (n±, . . . ,n* N ) 
with n\ = n* m = 1 is the unique partition of iV with ni > 0, n m > 0. Hence, the 
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measure p, such that p(rj*) = 1 can be viewed as a Gibbs distribution with arbitrary 
weights Cfyjv > 0, a m ^ > 0, all other weights equal to zero and the partition function 

Cn = (<fy,jV<2 TO| jv) 1 - 

In the rest of the paper we assume that the initial distribution is the total dissipation 
with a hN = ai jJV (0) = 1, A > 1, so that C w = A!, A > 1. 
Corollaries 

(A) The general form of weights dk,N(t)- 

The recurrence relation (12.41) allows to find a general form of the parameters dk,N{t), 
for an arbitrary function / > 0. To do this, we need two more notations. Let 

s=l 

be a partition from the set Qj of integer partitions of j and set 

s=l 

for j = 1,...,N and / > . 

Assertion. For given A and a function / > 0, the solution a/ Cj Ar(t) of (12.41) has the 
form 

a fclJV (t)= Yl Btf k hN)exp(-qtf k hN)i), t>0, k = 1 . . . , A, (2.15) 

where the coefficients B(r]( k >; A) do not depend on £ > 0. 

Proof is done by induction in 1 < k < A. Recalling our assumption that 



0>k;N '■— °fc;Af(0) 

we obtain from (12. 4p 



1, if k = 1 

0, if k = 2, . . . , A", 



a hN (t) = e-^-D/W*, t > 0, 

which is of the form (12TT5D with q(r]^;N) = (A - 1)/(1), B(r]^;N) = 1, A > 1. 
Consequently, we derive from (12.41) . 

2;7vl J 2(A — 1)/(1) — (A — 2)/(2) v /' " ' 
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which is of the form (I2.15P induced by the two partitions (2,0) and (0, 1) of 2. Let now 
(jZISD hold for 

a j,N, J = — 1. Substituting (I2.4p for a^jv^), j = 1, • • • , k — 1 into 

f)2.15p . the claim for a^A^t) follows from the fact that for any given 77^ G Qj : r]^ k ~^ G ilk-j, 

q(r] ij) ; N) + q(r] ik ~ j) ; N) = q(r] {k) ; N), 

where rj^ G ilk is the partition of k obtained by merging the two partitions r}^\r]^ k ~^. 

(B) Additive rates of coagulation: tf)(i,j) = (i + j)v, v > 0. In this case the 
probability distribution p(r], p; t) can be found explicitly. Namely, 

p( V , P ;t) = Me- N ^(l - e - N *) N - r N-( N ^ ]J t > 0, 

k=i 1 ' Hk ' 

{n u ...,n N )e fi r ,7v, (2.16) 

where Q r . t N = {v = ( n ii • • • 1 h n) G Qn '■ n i + • • ■ + n N = r , 1 < r < N} denotes the set 
of all partitions of iV with r summands (=clusters). 

Proof. The case considered conforms to ( 12. 3 p with f(j) = v > 0, j ' > 1. We will seek 
the solution of the recurrence relation (12.41) in the case considered, in the form 



a k , N (t) = e-( N - k ^ t {l-e- mt ) k - 1 v k , N , k = l,...N, N > 1, (2.17) 

with constants Vk,N > that will be determined from (12. 4p . For k = 1 (I2.17P reduces to 
(I2.12p with Vi } n = a ltN = 1. Next, by the induction argument, in k > 1 f)2.4p gives 

/ f* Q,kp~( 2N ~ k ) vu ('\ — P -Nvu\k-2 (N-k)vu N 

a) = / — ^ — ^ — d«) x; = 

i+j=k 







-{N-k)vt (-, -Nvt\k-i N l k ST^ 



2(fc 

In view of ( I2.17p . this leads to the following recurrence relation for Vj-^n' 

N^k 



VhN 



— 77 V hNVj,N, k = 2, . . . , N - 1. 



This is a particular case of the recursion (2.33) in [II], whose solution is 

7 & — 1 

v k ,N = N~^-—, k = l,...,N 
k\ 
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(see (2.39) in [II]). ■ 
From ( 12.16)) we find the probability p CO ag,Nif) of coagulation into one giant cluster of size 
TV by time t > 0: 

PcoagAt) = (1 - e"^)^ 1 -> 1, JV OO, t > 0. 

This says that a strong gelation phenomenon holds at any time t > 0, as N — > oo. 

We discuss below the probabilistic meaning of the distribution (I2.16p . We rewrite 
MP as 

prO^^^-^^n^S' (2 - 18) 

(ni, • • • , n„) G fi r ,Ar, 

where B"^ = N(r - 1)\(N - r)\N^ N - r \ It was noted in [H] that the sequence of 
weights ^rp, A; > 1 is a particular case of weight sequences satisfying the Gnedin-Pitman 
condition of exchangeability of Gibbs set partitions [9] . It is easy to see that in ( I2.18P the 
term 

1 V-iV(r-l)«t/ 1 _ e -Nvt\N-r 

r — 1 J ^ ' 

expresses the transition probabilities of a pure death process with rates fi rj N — (r — 
l)vN, r = 1, . . . , N, as in (2.17) of [H]. Note that the aforementioned transition prob- 
abilities are binomial distributions with time-dependent probabilities of success, while 
B r N in (12.181) is the (N,r) partial Bell polynomial (see (2.37) in [H] for more details). 
In the context of our model (Br^) -1 serves as the partition function for the conditional 
distribution F(X^\t) = rj I = r) which is microcanonical Gibbs distribution on 

f2 r Ar, with weights prp that are independent on N and t. In view of this, we conclude 
from ( I2.18P that the distribution of the number of clusters := ni(t) + . . . + riN(t) 

at time t is binomial: 

F(\X%\t)\ =r)= {^Z^) e ' N{r ' 1)Vt { 1 - e~ Nvt ) N ~ r > t ^ °- ( 2 - 19 ) 

This fact was originally proven by Lushnikov (see (49) in [18]). It is interesting that in the 
case of the Marcus- Lushnikov CP with the iV-dependent additive kernel, the distribution 
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of |xjy^(i)| is also binomial, but with a different parameter. The latter distribution was 
derived by Aldous in [1] , from the interpretation of the process as the vector of sizes of 
the continium random tree. 

(C) Multiplicative rates of coagulation: if)(i,j) = 2ij. Correspondingly, f(i) = 
i, i > 1, so that under the assumed initial distribution, ( 12. 4p becomes 



a hN (t) = e- iN - 1)l , t>0, N>1, 
a k , N (t) = e-^ kt f e k2u V (ia hN (u)e mu ) { 3 a hN {u)e N ^)d Ul k = 2, . . . , N, t>0. 

JO .■ , „_7„ 



i+j=k 

Denoting 

ia i>N (u)e Nm = b i>N (u) 
we arrive at the following recurrence relation: 

bi tN (t) = e\ t>0, N > 1, 

b ktN (t) = ke kH [ e~ k2u k iN (u)b jtN (u)du, k = 2, . . . , N, t>0. (2.20) 
Jo i+j=k 

The important fact is that the function bi^(t) and, consequently, by virtue of ( I2.20p . 
the functions bk^if) do not depend on iV: &fc,jv(i) = bk(t), t > 0. Thus, in the case 
considered, 

p( V , p- 1) = N\e~ NH f[ V = (ni, • • ■ , n N ) G Q N , t > 0, (2.21) 

fe=i nh ' 

where bk(t) are defined by (12.201) . The right hand side of (12.211) can be viewed as the 
Gibbs distribution on fi^r with the partition function N\e~ N 1 and the weig hts k > 
1, t > not depending on N. This allows to employ the known exponential relation 
(see for references e.g. [7]) between the generating functions H(t;x) and V(t;x) for the 



sequences hk(t) := ^r, k > 0, t > and := ^r^, > 1, t > respectively: 

= e y(t;:c) , V(t;z) = ^v k {t)x k , H{t-x) = ^^x k , t > 0. (2.22) 

fe>l fc>0 

We note that the radius of convergence of the power series H(t; x) is zero and that by 
virtue of the exponential relation ( I2.22p . the same is true for the series V(t; x). This says 
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that the two power series should be treated as formal ones. From the exponential relation 
(I2.22p it is easy to derive the following recurrence relation between the sequences {vk(t)} 
and {h k (t)}: 

n 

h {t) = l, (n+l)h r * 1 (t) = ^(k + l)v k+1 (t)h n - k (t), n = 0,l,..., t>0. (2.23) 

fc=0 

It goes without saying that f!2.23j) and (12.41) are equivalent. However (12.231) is much more 

convenient for the study of the asymptotics of v k (t), as h — > oo. Based on the fact that the 

k 2 t 

functions h k (t) = grow very rapidly with k for any fixed t > 0, we will demonstrate 
that the solution of (I2.23P is given by 

v k (t)~h k (t), k^oo, t>0. (2.24) 

Firstly we see from (12.231) that 

v n +i(t) < h n+1 (t), t>0, n = 0,l,..., (2.25) 

This implies 

n-1 

(n + l)v n+1 >(n + l)h n+1 (t) - J2( k + l)h k +i(t)h n - k (t), t > 0, 

fc=0 

and consequently, 



v n+ i(t) ^4 (k + l)hk+i(t)hn-k(t) 



> 



K+i(t) '" (n + l)h n+1 (t) 



This together with (12.251) say that for the proof of (12.241) one should validate the limit 

^ (k + l)h k+1 (t)h k (t) = ^ fn\ 2k{M)t n ^ 
(n + l)h n+1 (t) ^\k) 



k=0 y ' ^ y ' k=0 

for any t > 0. We have 

[n/2] / v [n/2] fc-1 



e~ 2n * 



J2[ n k ) e- 2fe(n - fc - x)t = e- 2 "* (l + ]T i - J>~ 2( "~ fc_1) *) < 

fc=0 ^ ' k=l ' j=0 

[n/2] 



-2«t + yne- 2 ^- 1 ^ -+ 0, n ->• oo. (2.26) 
fc=i 



15 



In view of the relation m = ( n ™ fc ) , the same limit can be proven for the sum 

Jt-l , V 



n-l 

-2nt 



e 

fc=[n/2]+l 

This completes the proof of (12.241) . Equipped with this result we are now in a position to 
study the time dynamics of clustering of groups of different sizes, as N — > oo and t > is 
fixed. With an obvious abuse of notation, let the random variable n kt N{t) be the number 
of groups of size k at time t > 0. By the known formulae for the functionals of the Gibbs 
distribution considered (see e.g. [7]), we obtain for any fixed moment of time t > 0, 

En ktN (t)=v k (t)^=^- t k = l,...,N. 
h N (t) 

cov{n KN {tln, N (t)) = v k (t) Vl (t) (^g^ - , k ^ I = 1,2 . . . , N. 

(2.27) 

Substituting h k (t) = it is easy to find that the three quantities tend to zero as 
— oo for any fixed k, I and t > 0. On the other hand, applying (I2.24p . we conclude 
that at any time t > 0, 

En N:N (t) -> 1, N ^ oo. 

This means that p CO ag,N{t) — > 1, as iV — >• oo at any moment £ > 0, which is equivalent to 
say that a strong form of gelation occurs during all time evolution of the process. 

Unlike this, the Marcus-Lushnikov process with iV-dependent multiplicative kernel 
exhibits gelation only after time t = 1. This fact was proven in [5]. Also note, that in [6] 
the aforementioned Marcus-Lushnikov process was represented as a random graph process, 
which among other things allowed to obtain a version of (12.211) for the case considered. The 
survey pQ enlightens this interesting connection to random graphs. Finally, we note that a 
properly time-space rescaled Marcus-Lushnikov CP with multiplicative kernel converges 
to a limit process called standard multiplicative coalescence (P-J). This fact facilitates 
the study of the emergence of the giant component. Regarding the formulae (12.271) . it is 
in order to observe that they are not valid in the general case of TV-dependent weights 
afe,jv(i). This can be seen from the derivation of (12.271) (see e.g. [7]). 
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(D) Becker-Doring pure coagulation process. This CP is a stochastic version of 
Becker-Doring kinetic equations proposed in 1935, to model a variety of phenomena in 
which only coagulations with monomers (=clusters of size 1) are allowed (see for references 
e.g. [3]). Formally 

fl>(i,j) = °> if min{2, j} > 1. 
Clearly, the rates of the process are of the form (12. 3p iff the function / has the form 

,0, if i > 1. 

/(<) = < (2-28) 
v > 0, if z = l, 



which leads to the rates 

V>(i,l) = V(M) 



u>, if z > 1 
2v, if z = l. 



Correspondingly, (12.41) takes the form 

a hN (t) = e -(tf-i>*, t > 0, 

a fclJV (t) = (k - l)v [ e^-^a^Mdu, 
Jo 

k = 2,...,N, t>0. (2.29) 

From (12.15P it is not difficult to derive that in the case considered a^^it) is a polynomial 
of degree k in z(t) = e - {N ~ l)vt , t > : 

k 

a k , N ^) = Y, m ^- e ~ i{N ' 1)Vt ^ k = l,...,N, t>0. (2.30) 

i=0 

By (I2.29p . the coefficients = m^N are defined recursively by 
(A; — l)TOj_i 

mi ifc = _ ' , i = l,..., k, k>2, m^ k = - 2_^m^ k , k>2, (2.31) 

with the initial conditions 

m ,i = 0, mi,i = 1, m , 2 = a 2 ,7v(oo) = 7^y^yj> m 2,2 = _ ^\TrT)- ( 2 - 32 ) 
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The recurrence relation in ( I2.3ip is validated by the induction argument: 

(k - l)v fe-^^(ym i|W • e^'^du = 
J ° V i=o 7 

(fc-l)t, T (Y, m ^- e ' l(N ~ 1)VU ) du = (^^'E n (l-e^- 1 )") = aUt), 

^° i=l i=l ' 

k = 2,...,N, t>0. 
Solving f)2.3ip . we find the explicit expressions for the coefficients m i fc : 

(-l)fc- 1 

m fe _i >fc = 0, fc > 1, m k ,k = _ k>l, 

Consequently, we see that the probability of appearance of one giant cluster of size N 
decays "almost" exponentially to zero, as N — > oo, at any time t > : 

PcoaMt) = N\a N , N (t) ~ iV!mo,iv = (i ^~^ 2 ~ V^(N - l)f e"^" 1 ), (2.34) 

where the first ~ is due to the last relation in (I2.3ip . 

In contrast to the models in (B) and (C), the Becker-Doring CP has a nontrivial 
equilibrium distribution (=measure) Hn(v)i which is given by the weights 0*^(00). We 
have 

k — 1 

a hN (oo) = 0, a k , N (oo) = m , k = _ ju=l » ^ > 2 > 

so that 

N h — 1 i 

Mv) = n\(n - i)-<*-w> Ui^Y—v ( 2 - 35 ) 

fe=i 

where we agree that 0° = 1 and \rj\ = ni + • • • + tin is the number of clusters in 77. It 
follows from (I2.35|) that [1n(v) = 0> f° r & H ^ = ■ ■ ■ 1 n N) G ^tv with ni > 0. 

Finally, we derive from (I2.35P the conditional distribution, given \q\, at the equilibrium 
of the process (=microcanonical distribution at equilibrium): 

N h — 1 1 

M , i w =i) n(^)"'^, 



k=l 
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r] = (m, . . . ,n N ) e Q N : = I < N, (2.36) 

where B^n = P(|t/| = I) is the (I, N) partial Bell polynomial on the set of all partitions of 
N with ni = 0, induced by the weights Vk ■= ^r- not depending on N. 

3 CP's with canonical and microcanonical Gibbs dis- 
tributions 

The distribution (12. 2p can be written as 

TV 

p(r), p; t) = C N (t) exp ( - ^(-Ttj log a ijN {t) + log(^iQ) j , ?7 e ^/v, * > 0. (3.37) 

i=i 

This shows that in the context of statistical physics, ( 12. 2 p conforms to the canonical Gibbs 
distribution with time dependent potentials of the following special form: 

TV 

H k ( V ;t) = 0, k>2 H 1 ( m t) = ^(-THloga i , N (t)+log(n i \)), 

i=i 

t>0, 7] e Q N . 

In the context of stochastic models of coagulation and fragmentation, Gibbs distributions 
(I2.2p with weights not depending on t and N emerged as early as in the 1970's in the works 
of Kelly (see [IS]) and Whittle (see [23]), devoted to reversible models of clustering at equi- 
librium. Vershik (see e.g. [23]) who intensively studied such Gibbs distributions in the con- 
text of equilibrium models of ideal gas, called them multiplicative measures. Pitman [21] 
introduced Gibbs processes of pure coagulation/fragmentatation on the state space fl[N], 
which is the set of all partitions of the set [1, . . . , N], while developing Kingman's theory 
of exchangeable partitons. Unlike the present paper, Pitman's definition of Gibbs process 
requires that microcanonical (rather than canonical) distributions are Gibbsian at any 
time t. Subsequently such processes were extensively studied by Berestecky and Pitman 
in [lj, the main result of which is the characterization of weights of Gibbs fragmentation 
processes on Q.\m. It turned out that the time-reversal of these processes are CP's with 
ifj(i, j) = a(i+j) + b. (Note that for b ^ the latter rates are not of the form (12. 2p ). In [TU] 
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Goldschmidt, Martin and Spano constructed a Gibbs fragmentation process with weights 
that does not obey the characterization condition in [3]. This became possible because 
the constructed process does not possess the mean-field property. The interplay between 
the set up, when the state space of a CP is set partitions, and the set up in the present 
paper is discussed in more details in [21] , [I] , [TT] . In [TTj . as a development of the idea 
of [TH], a characterization of coagulation- fragmentation processes, such that the induced 
birth and death processes are time homogeneous, was established. Based on this, 

a characterization of coagulation-fragmentation models, possessing time-independent mi- 
crocanonical Gibbs distributions ¥(X^\t) = n I = k) was obtained. By [11] and 

our theorem, the interrelation between CP's with Gibbsian canonical distributions and 
the ones with time-independent Gibssian microcanonical distributions is as follows. (We 
note that Gibbsian canonical distribution induces a Gibbsian microcanonical distribution, 
the latter being in general time- dependent). The CP's with ifi(i,j = a(i + j), a > are 
the only ones that have Gibbsian canonical and time-independent Gibbsian microcanical 
distributions. On the other hand, the CP's with ip(i,j) — a(i + j) + b, a > 0, b > 
are the only ones that have time-independent Gibbsian microcanonical distributions and 
non Gibbsian canonical distributions, while the CP's with = if(j) + jf{i), when 

/ ^ const are the only ones having Gibbsian canoninal and time-dependent Gibbsian 
microcanonical distributions. In the conclusion, we mention a representation of Gibbs 
distributions ( 12. 2 j) arising in the field of random combinatorial structures. In the case 
when the weights ai~,N(t) do not depend on t and N, the distributions ( 12. 2p depict the 
distributions of vectors rj = (m, . . . , Un) £ &n of component counts of random combina- 
torial structures (see e.g.[2],[8],[l],[T2],[T7]). In this set up nk stands for the number of non 
decomposable components (e.g. cycles in a random permutation) of size k. A cornerstone 
fact in the theory of random structures is the representation of the aforementioned mea- 
sures via the so called conditional relation, which proved to be very useful for problems 
of asymptotic enumeration. The version of the conditional relation for time dependent 
Gibbs distributions ( 12. 2 p reads as 

N 

p(r), p; t) = F[Z ltN (t) = ni, . . . , Z N)N (t) = n N I ^ kZ kjN (t) = N) 

k=l 
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rj = (n 1: . . . n N ) e f^Ar, t > 0, 

where {Zk^{t), k = 1, . . . , N, t > 0} is the triangular array of Poisson random variables 
with parameters cik,N(t), such that Zk t N(t), k — 1, 2, . . . , N are independent at any given 
time t > 0. We note that in our setting, 

N N 

F(j2 kz kAt) = N ) = {N^- l exp{ -J2a k , N (t)), t > 0. 

k=l k=l 
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